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An elastic string embedded in 3D space and subject to a short-range correlated random potential 
exhibits marginal pinning at high temperatures, with the pinning length L C (T) becoming expo- 
nentially sensitive to temperature. Using a functional renormalization group (FRG) approach we 
find L C (T) oc exp[(32/7r)(T/Td p ) 3 ], with Td p the depinning temperature. A slow decay of disorder 
correlations as it appears in the problem of flux line pinning in superconductors modifies this result, 
lnL e (T) ex T 3/2 . 



Elastic manifolds subject to a disorder potential 0] de- 
fine a rich and challenging problem in modern classical 
statistical mechanics, with numerous applications in con- 
densed matter physics jl],^ JjL[| . Fundamental questions 
arise regarding the relevance of the disorder [El and its 
quantitative effects on the manifold's quenched fluctua- 
tions on short || and large scales JjJ . Such static results 
then determine, via scaling arguments, the creep type 
motion of the elastic manifolds under the action of an 
external force, rendering these studies also relevant in 
the context of applied physics, e.g., plastic deformations 
in metals |3| or dissipation in superconductors Clas- 
sifying the problems through the dimensionality d of the 
manifold and the number n of transverse dimensions, the 
(1 Tridimensional problem describing strings moving in 
n directions, also known as the random polymer prob- 
lem, is particularly interesting |m. Here, depending on 
the number of transverse dimensions n, the temperature 
induced fluctuations compete in various ways with the 
fluctuations due to quenched disorder. In this Letter, we 
analyze the competition between elasticity, disorder, and 
temperature for the (1 + n)-string problem. In particu- 
lar, we present a functional renormalization group anal- 
ysis of the (1 + 2)-dimensional problem (a string in 3D 
space), where the disorder turns marginal at high tem- 
peratures; we determine the depinning temperature Td p 
above which thermal smoothing leads to a collapse of 
the effective disorder strength and the temperature de- 
pendent pinning length L C (T) beyond which the disorder 
dominates over the elasticity. 

The free energy F(u, L) = —TlnZ(\i,L) of a string 
starting and ending at the points (0, 0) and (u, L) [the 
first (second) coordinate denotes the displacement of (po- 
sition along) the string)] is given by the partition function 



(u,i) L 
Z = j V[u') exp {- I jf dz' 

(0,0) 



du' 



+ U(u',z') 



with e the elasticity of the string. The disorder potential 
U(u, z) is assumed to be a random gaussian variable with 
zero mean and a correlator 



with Kq(u) a smooth function decaying on the scale £. 
The partition function Z describes the competition be- 
tween the elastic energy (e/2)(<9 z ii) 2 > and the disorder 
potential U(u, z), from which the string can gain energy 
by choosing minima with U (u, z) < 0. 




(U(u, z)U(u', z')) = K a (\u - u'\)6(z - z'l 



(1) 
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FIG. I. Schematic picture of the FRG flow for the 
marginal case with n = 2: At temperatures below depin- 
ning (T < T dp w (e_Ko(0)£ 2 ) 1/3 ) the flow proceeds through 
the Larkin regime (with a fixed point characterized by the 
wandering exponent £ L = 3/2) towards the random mani- 
fold fixed point with its non-trivial exponent £2 ~ 5/8. 
At high temperatures T > Td p the flow starts out in the 
thermal region with fth = 1/2. The linear flow towards 
the 'Boltzmann' fixed point (BFP) rapidly erases the details 
in the initial condition and produces a gaussian correlator 
Ki(u) ~ gi exp(— u 2 e/47r£T) with weight gi. The non-linear 
terms in the FRG produce an overall growth of the correla- 
tor while the linear terms keep its shape to the 'Boltzmann' 
form (adiabatic pumping). The sharp increase of the prefac- 
tor gi at l c defines the pinning length L C (T), beyond which the 
flow enters the random manifold regime with a rapid change 
of the thermal exponent to the random manifold exponent 
(£th = 1/2 — > C(0 ~~ > C2 ~ 5/8) and a renormalization of the 
temperature towards zero. The flow ends in the zero temper- 
ature random manifold fixed point. 

The quantity characterizing the behavior of the string 
is the displacement correlator (u 2 (L)) = ([u(L) — u(0)] 2 ) 
oc L 2 ^ describing its wandering with distance L. We dis- 
tinguish between the perturbative (Larkin j^]) regime at 
short scales L < L c and the random manifold regime 
H at L > L c , as well as temperatures below and above 
the depinning temperature Td p ~ (ei4T (0)^ 2 ) 1 / 3 , see Fig. 
[l]: For T < Tdp the wandering exponent £ L = 3/2 at 
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small distances L < L c (see [||), while at large distances 
L > L c , depends on the dimensionality of the space: 
Ci = 2/3 is an exact result [Q while for n = 2 numeri- 
cal simulations || give a value C2 ~ 5/8. The crossover 
length L c separating these two regimes has the form Q 
L c w [e 2 C 4 /i^o(0)] 1/3 . At high temperatures T > T dp 
the short scale behavior is subject to thermal smoothing, 
resulting in a thermal wandering with £ t h = 1/2, while 
at large scales L > L C {T) the disorder induced line wan- 
dering characterized by £ n prevails. The crossover length 
L C (T) increases algebraically with temperature L C (T) ~ 
L c (T/T dp ) 5 for the (1 + l)-case §, while for the (1 + 2)- 
dimensional problem L C (T) is exponentially sensitive to 
temperature with @ L C (T) ~ A(T) exp[C(T/T dp ) 3 ], 
and C a constant of order unity. It turns out that the 
usual perturbative methods neither provide the value of 
C nor that of the prefactor A(T). The exponential sen- 
sitivity of the pinning length L C (T) to the temperature 
T is related to the appearance of a phase transition in 
dimensions n > 2, separating a low temperature disorder 
dominated phase from a high temperature thermal phase 
characterized by (th PC - The (1 + 2)-problem then cor- 
responds to the lower critical dimension for this phase 
transition and thus exhibits marginal behavior. 

The main goal of the present work is to develop a con- 
sistent scheme for calculating the pinning length L C (T) 
in the (1 + 2)-problem specifying the coefficient C in the 
exponent. While this goal can be achieved within a one- 
loop calculation including three-replica terms, it turns 
out that the determination of the prefactor A(T) requires 
a two-loop analysis (accounting for four-replica terms) 
which is beyond the scope of this letter. We first discuss 
perturbative methods for calculating L C (T) at high tem- 
peratures and different dimensions n. These techniques 
fail for the (1 + 2)-problem and we will then turn to the 
more powerful renormalization group method. 

Perturbation Theory: With the weak disorder pro- 
viding the small parameter A = J d n u Kq(u) we ex- 
pand the fc-fold replicated Green function (see Jl2| ) 
G(ui, . . . , us,-; L) = (Z(ux,L) . . . Z(uk,L)) and calcu- 
late the mean squared displacement (u 2 (L)) = limfc_>o 
/ d n ud n U2 ■ ■ ■ d n Uk u 2 G(u, u 2 . . . , Ufe; L). The disorder is 
relevant when (u 2 (L — > 00)) grows faster than the ther- 
mal wandering (u 2 (L)) t h = (nT/e)L which is dominant 
at small lengths; the comparison of the two expressions 
then defines the characteristic crossover length L C (T), 
(nT/e) L C (T) ~ (u 2 (L c (T))). To lowest order the mean 
squared displacement takes the form 



(u 2 (L)) 



nTL 



(A/T 2 )( £ /T)»/ 2 W2 
(47r)™/2(4_ n ) 



(2) 



for n < 2 the disorder corrections grow more rapidly 
than the thermal wandering and we obtain Q L C (T) ~ 
L c [T/(eA£) 1/3 ] 5 for n = 1. For 2 < n < 4 the contri- 
bution from the random potential goes to zero at large 



scales and the disorder is irrelevant |yj . The case n = 2 
is marginal: to lowest order the disorder correction is L- 
independent and the next order term oc A 2 provides a 
logarithmic correction, 



(u 2 (L)) = 



2TL 



1 eA 



1 



2 A 2 



16tt 2 T 6 



ML/0 



(3) 



producing an exponential temperature dependence 
L C (T) oc exp[CT 3 /eA] with C of order unity (we com- 
pare terms oc A 2 and oc A 3 ). The arbitrariness of the 
criterion prevents us from finding the coefficient in the 
exponent and the prefactor and we have to develop a 
more systematic way in order to deal with the problem. 

Renormalization group: The basic idea is to construct 
the renormalized effective correlator Ki of the disorder 
potential describing the behavior of the manifold on short 
and large scales. With the disorder becoming irrelevant 
for (internal) dimensionalities d > 4 || , the effect of dis- 
order is analyzed within an e = 4 — c?-expansion. Start- 
ing from short scales, the RG flow goes through a special 
point where Ki becomes singular. This point is iden- 
tified with the pinning length L c of the string p3[ . Its 
physical meaning is clear: the fact that the correlator be- 
comes singular at L c implies that the perturbation the- 
ory breaks down and the manifold splits up into indepen- 
dently pinned domains of size L c - By way of introduction 
we briefly derive the zero-temperature pinning length L c , 
starting from the one-loop FRG equation HQ 

diK t = (4 - d - 4C - n()K t + CV • (uK t ) 

+ I[K^ v {u)K^{u)/2 - K^(u)K^(0)}, (4) 

with £(/) = (1/2)3 \n(u 2 )/d In L the wandering exponent, 
/ = A d /(2ir) d e 2 A i - d {A d = 2ir d / 2 /T(d/2), and A" 1 - £ 
is the short scale cutoff), and the indices /z and v denote 
derivatives with respect to the Cartesian coordinates 
and it„. The RG variable I is related to the length L via 
/ = In AL. We proceed with differentiating ([|) four times 
with respect to u K . Substituting u = 0, the wandering 
exponent £ drops out and we obtain 



a i r ; = (4-d)r ; + [/(n + 8)/3]r 2 , 



(5) 



where T; = d 4 Ki/du^\ u= Q. Integrating (||) we encounter 
a divergence at l c ~ ln[3(4 - d)/I(n + 8)r ]/(4 - d) 
(we assume weak disorder with LTq <C 1), allowing us 
to define the collective pinning or Larkin length L c — 
A" 1 exp(; c ) « [e 2 ^ 4 /if(Q)] 1 /(4-<i) w here simple pertur- 
bation theory breaks down. The FRG equation (|J) al- 
lows us to further characterize the Larkin regime: On 
short scales we can neglect the non-linear term and an 
expansion of the correlator Ki(u) ~ a(l)u 2 /2 produces 
the wandering exponent £ L = (4 — d)/2 at the Larkin 
fixed point Note that, although Eq. (g) is written 

in the form of an e-expansion, the result for L c is valid 
for any e as long as we investigate the short scale behav- 
ior, where the functional renormalization group is just 
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a way to sum the perturbation expansion. Indeed, one 
can show that in the non-marginal situation higher order 
loop corrections to Eq. ^ do not change the result for 
L c , while in the marginal case (e.g., for d = 4 at T = 0) 
we have to account for two-loop corrections. 

Next, we proceed to study finite temperatures, 
where besides a simple diffusive term oc TAKi, new 
three- and higher replica terms are generated describ- 
ing non-gaussian fluctuations in the disorder potential 
U. Generalizing the disorder statistics to include non- 
gaussian terms and replicating n times we arrive at 
a new Hamiltonian TL n — 7if rcc + 'Hint with 7Yf rec = 
I d d zJ2 a (e/2)(d z u a ) 2 and 



Hi, 



d d 2 



a, 13 



2T 2 



+ . 



Here, u a/ 3 = u Q - U(j and S(u, u) is the three-replica 
term describing non-gaussian fluctuations in the disorder 
potential (the . . . stand for higher-replica terms depend- 
ing on three or more variables). At finite temperatures, a 
consistent analysis to one loop order requires to include 
both two- and three-replica terms K (u) and S(u, u) in 
the FRG flow. After a second order cumulant expansion 
in the two- and three-body interaction Ti.- ln t, we proceed 
with the standard momentum shell RG, expand the u- 
fields into Fourier series and perform the integration over 
fast modes. The resulting equations take the form 

diK,(u) = (4-d-4C-nO-Ki + CV • (uK t ) + I'TK^ 

+ w/2)Kr ( U )/2 - xr^Krm 

- 2I'T[S^(0, u)+Sf*(u, 0)-[9 u -a a ^](u, u) 

+ [9 u -a a 5 / ](o,u) + [a u -a a 5 / ](u,o)], (6) 

diSi(u,u) = (6-2d-6C-2nC)5, + C[V-(u5,) + V.(u5i)] 

(7) 

d l T l = (2-d-2QT l , (8) 



I'T[S^ 
ITK^{vl)K^{vl)/A, 



where I' = IeA 2 and the superscripts /i (/2) denote deriva- 
tives with respect to the first (second) variable in S/(u, u) 
(implicit summation over double indices is assumed); 
similarly, the gradients V (V) denote derivatives with 
respect to u (u). The system (||) - (||) has to be solved 
with initial conditions Ki = q(u) = Kq(u), <5';=o(u, u) = 
(the bare potential is gaussian), and Ti = o = T. Note that 
the three-replica term is driven by a source oc TKf and 
hence is irrelevant in the T — analysis. Furthermore, 
four- and higher replica terms can be omitted within 
the present one-loop analysis: with Ki(u) = 0(A) the 
three-replica term S/(u, u) = 0(A 2 ), see (R), while the 
four-replica term is driven by a term oc Kf and hence 
-D(ui,u 2 ,u 3 ) = 0(A 3 ). However, going to second-loop 
order the four-replica term Di has to be included. 

The above system (|J) - (|^) can be solved asymptoti- 
cally exactly for weak disorder and we proceed with the 



analysis for the elastic string, d = 1. On short scales the 
behavior of the string is dominated by thermal fluctua- 
tions and the wandering exponent is equal to £th = 1/2; 
as a result, the temperature is not renormalized, see 
Eq. (||) . On the other hand, on large scales C approaches 
the value (2 ~ 5/8 and the behavior of the string is gov- 
erned by the zero temperature fixed point. It turns out 
that the temperature starts renormalizing only close to 
L C (T), allowing us to fix £ = 1/2 in Eqs. (g) - (§). The 
non-linear and three-replica terms are of second order 
in the disorder strength and thus unimportant at small 
scales I < ! c , allowing us to neglect the terms oc Kf and 
oc Si in (||) in the initial stage of the flow. The Ansatz 



Xi(u)=i^(u)exp[(l-n/2)/], 



(9) 



then maps (g) to the Fokker-Planck equation describ- 
ing the 'time' evolution of the probability distribution 
function Pi for an overdamped particle in the harmonic 
potential V(u) = u 2 /A at a 'temperature' T' = T/ireA : 
diPi = V • (uP/)/2 + T'APi. At high temperatures T, 
the distribution function Pi(u) rapidly approaches equi- 
librium and takes the form of a Boltzmann distribution, 



P, 



A(47rTT™ /2 exp(-u 2 /4T'). 



(10) 



The function P B {u) depends on the disorder potential 
only via the integral weight A = J d n u Kq(u), i.e., the 
evolution erases the details of the initial condition. De- 
pending on the dimensionality n, the following situations 
arise (c.f., (||)): For n > 2 the solution of (||) decreases 
with increasing I and the disorder renormalizes to zero, 
in agreement with the existence of a high temperature 
phase when n > 2, see pT[ . For n < 2 the correlator 
Ki(u) = exp[(l — ?i/2)l]P B (u) grows (until reaching the 
size of the non- linear term) and the disorder is relevant. 
Comparing terms linear and quadratic in Ki we find the 
collective pinning length L C (T) ~ L c [T/(eA£) 1/3 ] 5 . Most 
interesting is the marginal case n — 2 where we find 
a 'Boltzmann' fixed point. In order to obtain L C (T) 
we need to take into account the terms 0(A 2 ) effec- 
tively pumping additional 'particles' into the system and 
leading to an overall growth in the disorder strength 
A/ = J d n uKi(u). However, for small A; the strong dif- 
fusion will relax this 'pumping' towards the 'Boltzmann' 
distribution, i.e., the solution of Eq. (o) keeps the form 



Ki(u)= 9l (T /I) exp(— u 2 /4T'), 



(11) 



with the dimensionless coupling (7; = (7r/4)A/e/T 3 . After 



diagonalizing the term [S, 



d u -duSi] in (0) we 



find its Green function G s (u, u; u', u'; I) = l/12[nI'T(l- 
e- 1 )} 2 exp{-[(u-u'e- i / 2 ) 2 +(u-u'e^'/ 2 ) 2 +(u-u' e -'/ 2 ) • 
(u-u'e~' /2 )]/3/'T(l -e~')} and obtain the three-replica 
term Si (u, u) by simple integration. Substituting Si back 
into (^) and integrating over u we find that the three- 
replica terms mutually cancel and thus do not contribute 
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to the flow of gi. The amplitude gi then is determined 
by the remaining non-linear terms alone, 

dm = gf/8, (12) 

with the solution gi = go /{I — gol/8) as derived from the 
initial condition g = (7r/4)(eA/T 3 ). At l c = 8/g Q the 
amplitude gi c w 1 and we find the pinning length 

L c (T)~exp[(32/7r)(T 3 /T 3 p )], T 3 p = e A, (13) 

the main result of the paper |ll| : the one-loop FRG anal- 
ysis allows us to fix the coefficient C in front of the pa- 
rameter l/go in the exponent. Although the criterion 
gi c ~ 1 might seem arbitrary it is sufficient to find the 
correct asymptotic value of the RG variable l c as long 
as the initial value go is small. Indeed, the solution gi 
remains small in most of the interval [0, l c ] and grows 
rapidly only very close to l c . 

An interesting alternative is offered by mapping the 
random polymer problem to the Kardar-Parisi-Zhang 
equatio n p7[ , where the dynamic renormalization group 
(DRG) ]18| provides excellent results for the n = 1 prob- 
lem jL9| . The weight g in the (white noise) Ansatz K wn = 
gS 2 (u) replacing ( p"l| ) obeys the equation Q ®vg = (2 — 
n)g+[2(2n-3)/n\§*, with g(0) = (A n /(27r) n )eA/T 3 (the 
prime indicates scaling along the transverse dimension). 
For n = 2 the solution 'explodes' at l' c = 27r(T/T dp ) 3 , 
implying that L C {T) ~ exp[47r(77T dp ) 3 ] ||, different 
from Eq. @. Note that, in contrast to the FRG, the 
DRG scheme does not renormalize the correlator func- 
tion. Here, we argue in favor of our result ((l^): While the 
DRG technique produces accurate results (including a 
random manifold fixed point) in dimensions n < 3/2, the 
results are less convincing in larger dimensions n > 3/2, 
where the fixed point gives way to a divergence. On the 
other hand, the FRG technique features both an appar- 
ent divergence at l c and a random manifold fixed point 
as I — ► oo in all dimensions n < 2, allowing for a consis- 
tent physical interpretation of the results in terms of the 
diagram sketched in Fig. 1. 

Above we have assumed that the integral J d n uKo(u) 
converges, however, flux lines in superconductors exhibit 
a long range interaction with the disorder potential spoil- 
ing this assumption: for values u exceeding the coherence 
length £, the correlator obeys the asymptotics K (u) ~ 
K a (0)(e/u 2 )ln(u/O and the integral J u<u , d n u K (u) 
diverges as K (0)£ 2 ln 2 (w*/£)- We estimate L C (T) for 
this case: Neglecting the non-linear term in the first stage 
of the RG transformation, we approximate the solution 
of the Fokkcr-Planck equation in the region u 2 < T£/e 
by [A ; /47rT']exp[-w 2 /4T'] [c.f., Eq. ©], with an l- 
dependent function A;, a consequence of the divergence 
of the integral J d n u Kq(u). The weight A; collects those 
particles reaching the central region on scale I: a displace- 
ment u* ^> \/T£/e of the harmonic oscillator takes a 
'time' 1* ~ ln(u*/^/T£/e) to reach the region u < y/T£/e 



and hence A/» ~ Ko(0)£ 2 l ■ The dimensionless cou- 
pling gi then grows already in the linear regime due to 
the 'pumping' from large distances. The pinning length 
follows from the condition g; c = eA; c /T 3 w 1 and we 
obtain L C {T) ~ exp[C(T 3 /eX(0K 2 ) 1/2 ] with C of order 
unity. The same result follows from a perturbative treat- 
ment. 

In order to find the prefactor A(T) in the expres- 
sion ( |l3| ) for L C (T) we have to account for 0(A 3 ) cor- 
rections, resulting in a flow digi — gf/8 + ngf with 
n > a constant. This equation explodes at l c = 
8/go — 64?7ln g , producing a prefactor of the form A(T) ~ 
A _1 (T 1 dp/7 1 ) 192,? . While the two-loop contribution n 2 -\ = 
275/5184 to n can be found from jl3|], the part arising 
from the four-replica term is more difficult to obtain — 
whether the four-replica terms mutually cancel (e.g., due 
to some symmetry) , as was the case in the one-loop anal- 
ysis above, remains an interesting and open problem. 
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